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There has been much recent study on the application of spin chains to quantum state transfer and communi- 
cation. Here we discuss the utilisation of spin chains (set up for perfect quantum state transfer) for the knitting 
of distributed cluster state structures, between spin qubits repeatedly injected and extracted at the ends of the 
chain. The cluster states emerge from the natural evolution of the system across different excitation number 
sectors. We discuss the decohering effects of errors in the injection and extraction process as well as the effects 
of fabrication and random errors. 
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I. INTRODUCTION 

With conventional information processing and communica- 
tions, optics provides the bandwidth and robustness for long 
distance communication. However, for communication over 
short distances, for example within and between adjacent sil- 
icon chips, signals remain electrical, to avoid the energy and 
cost overhead of conversion of information between different 
physical embodiments. Similar thinking exists in the quantum 
arena. Whilst quantum states of light are widely regarded as 
the vehicle of choice for quantum communication over large 
distances, there has been much recent interest in the poten- 
tial use of spin chains for quantum communication over much 
shorter distances. When the task at hand is communication 
within a quantum processor, or communication between ad- 
jacent processors or registers, it may well be that a chain of 
spins — the same hardware from which the processors and reg- 
isters are constructed — can play an effective and useful role 

ma. 

In its simplest guise the term "spin chain" applies to any set 
of two-state quantum systems coupled to their nearest neigh- 
bours. Clearly qudits or even continuous variable oscillators 
could replace the qubits, but as most quantum information 
studies generally focus on qubits, most spin chain studies do 
likewise. A chain could literally comprise spins or magnetic 
moments, such as with a string of fullerenesJ3|] or magnetic 
particles |01 or nuclear spins in a molecule H7 But it could 
also describe a system of electrons or excitons |6|-[8l in a chain 
of interacting quantum dots, or other devices. 

If the ground or prepared state of a spin-( 1 /2) chain is all 
spins down (|0)), then a complete single-qubit excitation (| 1)) 
is made by flipping one spin up. An arbitrary qubit state can 
thus be injected into a spin chain by preparing the "injection 
site" qubit in the appropriate superposition of up and down. 
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For quantum communication, the question is as to how well 
this state transfers — in terms of the fidelity of the initial state 
against that which emerges at the "extraction site" at a later 
time, resulting from the dynamics of the chain. Usually the 
injection and extraction sites are the opposing ends of a chain. 
Study of state transfer has been made for unmodulated chains 
[Toll , systems with unequal coup lings ifTlll . systems with 
controlled coupling at the ends Il2ll and parallel chains lf]~3tl . 
Of specific interest to our work here is the case of linear chains 
where the nearest neighbour couplings J^j+i between spin 
sites i and i + 1 are engineered to effect perfect state transfer 
(PST) between the injection (i = 1) and extraction (i = N) 
sites JH [Till. For a chain of N spins the PST couplings are 
given by UM 

J M+ i = J oy /i(N-i) (1) 

where Jo is a coupling that characterises the whole system and 
sets the timescale for PST, or mirroring, as t^j = Trh/2 Jq. 

In our work here we utilise such PST spin chains for a dif- 
ferent purpose — the construction of distributed cluster state 
structures. One potential application of short range quantum 
communication is to build up distributed entangled resources, 
that can be used, or consumed, to subsequently enable dis- 
tributed quantum processing through the concept of one-way 
computation lfl5ll . In this approach, the cluster state entangled 
resource lfl6ll is consumed by a sequence of measurements to 
effect the computation. Here we focus on the construction of 
entangled resources, which basically emerge from a suitable 
qubit injection and extraction protocol and the ability of PST 
spin chains to produce two-qubit entang ling gat es when oper- 
ated across different excitation sectors IU7l - ll9ll . For such an 
application, there is clearly merit in being able to generate the 
entangled resource as rapidly and effectively as possible. We 
shall demonstrate that an element of a cluster state ladder can 
be knit in time t^j, independent of the length N of the chain 
ll20ll by suitable injection and extraction of qubits at the ends 
of the chain. The construction of cluster state resources using 
PST spin chains as an entangling bus, with assumed access to 
all spins in the chain, has been discussed in lfl8l fl9ll . Here 
we adopt the original concept of a spin chain, where access is 
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restricted to the ends of the chain, and build the resource with 
a suitable injection and extraction protocol. We also consider 
the effects of realistic forms of decoherence acting in the spin 
chain and errors in the injection and extraction protocol. 



II. SPIN CHAIN DYNAMICS 

We start with the time-independent Hamiltonian that de- 
scribes the natural dynamics of a length JV nearest-neighbour- 
coupled spin chain 



JV 



N-l 



H = Y,E i \l)(l\ i + J2 Ji,i+i[\l)(0\i®\0)(l\ i+ i + 

i=l 

|0><1|<® |l)(0| m ]. (2) 
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We assume that the single site excitation (from down, |0), to 
up, 1 1)) energies Ei are independent of the site i, or are tuned 
to be so. The couplings J^+i are given by (|T). Tuning of 
the energies and couplings, such as via local fields or man- 
ufacturing control, is required, as per the PST scenario ((T). 
The total number of excitations T in the chain is given by the 
expectation value of the operator 



JV 



T = ]T|l)<l| 



(3) 



We also define the mirror operator M as that which reflects the 
state of the chain about its midpoint (which is spin (JV + l)/2 
for odd JV and the "gap" between spins N/2 and (N/2) + 1 
for even JV). So operationally M effects the following to each 
term in any arbitrary superposition state of the chain: 

M|o>i|6>a...|y>JV-i|«>JV = \z)i\y) 2 ...\b) N -i\a) N . (4) 

Clearly, for parameters restricted to achieve PST, the 
Hamiltonian (fJJ commutes with both T and M and so the 
system energy eigenstates \£f/) are also eigenstates of both T 
and M. For a chain of size N there are 2 N eigenstates in total, 
with T ranging from zero to N and each sector T containing 
JV!/(JV - T)\T\ eigenstates. It is also helpful to define the 
total chain "spin flip" operator T by: 



N 



^ = n (ii><oi, + io><iio 



(5) 



The eigenstates for excitation number N — T then follow from 
those for T by application of T . Within this framework it 
is straightforward to understand how PST, or more generally 
state mirroring rf2~ll I22I . occurs. Any initial state of a spin 
chain 1^(0)) can be decomposed into its even and odd (under 
M) parts, so 



|*(0)) = -^(|* + (0)> + |*-(0)>) 



(6) 



with |*±(0)) = ^(|*(0))±M|*(0)». Clearly the M 
eigenstates can be decomposed as superpositions of even and 



odd energy eigenstates ^±(0)) = Yl+k c ±fcl e ±fc) and then 
for the evolved state at time tu to have unit fidelity against 
the mirrored initial state M|^(0)), it must be of the form 
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It is therefore clear ifell |22| that quantum state mirroring 
places a requirement on the chain energy level spectrum so 
that the phases in the evolved state conspire to give the form 
© at the mirror time tM, with the coupling choices given in 
([D forming an example lfl4tl that produces a suitable energy 
level spectrum. 

The overall phase 9 in the evolved state (0 is potentially 
a hindrance when it comes to PST. In addition to its time 
dependence, this phase depends on both the chain length 
TV and the excitation number T. For the case of a single 
excitation (T = 1) the phase factor at the mirror time is 
exp(— i9(tM)) = (—i) 1 *^ 1 ITill and so it is recognised that 
(JV — 1) needs to be a multiple of 4 to eliminate it, although 
clearly in practice as JV will be known for a system where the 
couplings have been engineered to, for example, satisfy (fl3, 
this phase is a known correction, rather than unknown deco- 
herence. 



m. TWO-QUBIT GATES 

When it comes to remote quantum gates, which form the 
basis for construction of distributed entangled resources, this 
phase is the enabling effect lfl7l - [l9tl . For a given spin chain 
of length JV, PST does not in general work for superposi- 
tions that include states from different T sectors, due to the 
T-dependence of 9, although it works for arbitrary super- 
positions within a fixed T sector. However, this effect can 
be turned to advantage. Consider PST spin chains where 
qubit states are injected onto the two extremal spins of the 
chain (i = 1 and i = JV) at time t = 0. The initial 
state is thus in general a superposition of the T = 0, 1 
and 2 sectors. For the simplest case of the diagonal ener- 
gies in d2J equal to zero (Ei = 0) and expressed in the basis 
{|0)i|0)jv, |0)i|1)jv, |l)i|0)jv, |l)i|lM, at time t = t M the 
natural evolution of the chain effects a gate G given by 
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on the initial two-qubit state. We note that for the specific case 
where (JV — 1) is a multiple of 4 the T = 1 sector phase fac- 
tors are unity ltl4ll . but for all values of JV the gate G is a max- 
imally entangling gate, producing a concurrence of unity from 
an initial two-qubit product state given by an equal weight su- 
perposition of all four basis states. Thus the natural evolution 
of the PST spin chain enables a remote maxim ally -entangling 
two-qubit gate between the ends of the chain II17H19I1 . which 
cycles between four different variants with increasing JV. Un- 
derlying all these gates is an effective phase flip that ensues in 
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FIG. 1: EoF vs. rescaled time t/tM for spin chains of length TV 



FIG. 2: Entropy S vs. rescaled time t/tM for spin chains of length 
N = 9, 13, 17. 



the doubly-excited sector, which can be understood as arising 
from the anticommutation of two non-interacting fermions as 
they pass through each other IU8II . This picture is helpful in 
the multiple-excitation sector, where a phase factor of e~' lv 
results from every fermionic crossing. Experimentally, this 
can also be achieved via crossing beams of Rydberg atoms 
meeting in cavities l23ll . which therefore represent a potential 
alternative medium for this procedure. 

The dynamical two-qubit gate underpins the construction 
of the cluster state resource, so we first examine this basic 
gate, acting upon qubits each prepared in the state |+) = 
-^j(|0) + |1)) and injected simultaneously at the ends of the 
chain. The entangling action of the two-qubit gate is illus- 
trated in Fig. Q] where the time evolution of the Entanglement 
of Formation (EoF) ll24ll of the two end qubits is shown under 
the action of the PST spin chain dynamics. The EoF reaches 
unity at t,M = irh/2Jo and every irh/ Jq thereafter. We see 
that with an increasing number of spins in the chain, the width 
of the EoF peak decreases. However this decrease is not lin- 
ear, instead the dependence of the width of peaks in Fig. [T] 
with TV is well approximated by 1/ y/N within the range of 
chain lengths explored, demonstrating that even for very long 
chains recovery of the entangled qubits should be feasible. 

A single entangling gate essentially works for any length TV 
of a PST chain. However, in order to utilise additional qubit 
injections at the end of the chain, the length needs to be such 
that subsequent injections and extractions can be made inde- 
pendently. In order to understand this, it is helpful to consider 
the entropy of the end two qubits as a function of time. In 
order to calculate the EoF of the end two qubits, the spins 
comprising the rest of the chain (i = 2 to TV — 1) are traced 
out, to leave the density matrix p\ ,jy of the two end spins. If 
this is not pure, it provides a signature of entanglement be- 
tween the two end spins and the rest of the chain. This is 
illustrated in Fig. [2] which shows the evolution of the entropy 
S = —Ti^i^n log 2 Pi.n) as a function of time. 

Clearly the entropy is zero at the initial injection time and 
all integer multiples of tM- Of interest to us here is what hap- 



pens in between. It can be seen from Fig. |2]that provided that 
TV is sufficiently large, the end two qubits disentangle from the 
rest of the chain, whilst the excitations are propagating and lo- 
calised entirely in the middle region of the chain. The latter 
can be appreciated by looking at the site occupation probabil- 
ities, see Fig. [9f a). This results in S — > for windows (that 
widen for increasing TV) centred on odd half integer multiples 
of tM- Similar to Fig. [T] the width of the dip at t = tM 
decreases approximately as 1/ VTV within the range of chain 
lengths explored. For chains of length TV > 9, it is possible to 
independently inject further qubits. This will form the basis 
of the cluster ladder construction. Before considering the in- 
jection and extraction protocol, though, we first consider the 
effects of various errors on the basic entangling gate and the 
evolving two-end-qubit entropy. 



IV. ERRORS AND DECOHERENCE 

For non-zero but site-independent excitation energies for 
the qubits in the chain (the diagonal terms in Eq. (O), ad- 
ditional phases arise in the chain dynamics, but these do not 
affect the entangling capacity of the remote two-qubit gate. 
Of more interest is the effect of variations, or errors in the site 
energies. A simple estimate of the effect of errors in the en- 
ergy level spectrum ITil suggests an overall error, or loss in fi- 
delity for PST that scales linearly with TV. Figure[3]illustrates 
the EoF evolution as a function of the spin chain length TV, in- 
cluding a random modulation of the on-site energies Ei. The 
on-site energies are now given by Ei = eri, where < < 1 
is a random number from a uniform distribution. Each point 
in Fig. [3]corresponds to the EoF at t = tM averaged over 100 
random realizations. Figure [3] confirms that the loss in EoF 
scales linearly with the number of spins in the chain. More- 
over, it shows that the EoF decay is not linear for increasing 
s/Jmax, with J max = maxi{J itl+1 } = 1. For small and 
medium perturbations of the on-site energies, EoF close to 
unity can still be achieved for all chain lengths considered. 
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FIG. 3: EoF at t = tu vs. N for three values of e, as labelled. Data 
points are averaged over 100 random on-site energy realisations. 



FIG. 4: EoF at time t = tu for a 9-spin chain vs 7 and e. Data 
points are averaged over 100 random on-site energy realisations. 



As we are dealing with multiple excitations in the chain, we 
also consider the effect of interaction between excitations in 
nearby sites by adding the following term to Eq. (|2): 



N-l 

H'= ^ 7 J |l)(l| i ®|l)(l| i+1 . 

i=l 



(9) 



For example, H! may correspond to a biexcitonic interac- 
tion in quantum dot-based chains l25l l26ll . Figure [4] shows 
the dependence of the EoF at t = tu on the magnitude of 
the maximum on-site energy perturbation e and of the site- 
interaction perturbation 7 for TV = 9. Even with perturbations 
as big as 20% of the characteristic coupling strength Jo, the 
EoF hardly suffers and remains above 90%, indicating that the 
system is extremely robust against these sorts of errors. The 
influence of 7 in Fig. [4] is limited, as three quarters of the 
initial input state contain no more than one excitation and are 
thus not subject to the effects of interactions between excita- 
tions. 

Next we explore the effect on the remote gate of unwanted 
longer range interactions ll27i HH, which could be an issue 
when considering pseudospins based on charge degrees of 
freedom. We add to Eq. (O the perturbative term 
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with J M+2 = A(J M+ i + J i+ i )i+2 )/2 to simulate the origi- 
nal coupling modulation. A discussion about the experimental 
relevance of this term has been given in J29J on the example 
of graphene and self-assembled quantum dots. H" commutes 
with both T and M, so that a set of common even and odd 
eigenstates still exists. The dependence of the EoF at t = tu 
as a function of N is presented in Fig. [5] for three different 
values of A. The EoF displays a decay linear in N. For small 
values of A, relevant e.g. to graphene quantum dots ll29ll . 
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FIG. 5: EoF at t = tu vs. N for three values of A, as labelled. 
Lines are best fit to numerical data. 



the EoF is well conserved, although the effect is much more 
pronounced in very long chains as the increase in number of 
spins leads to more perturbation terms in Eq. ([ToV For exam- 
ple, with a value of A as large as 0.1, a short chain (N = 9) 
achieves an EoF of 0.79, with this dropping to 0.59 for a long 
chain (N = 30), as shown in Fig. [5] We also note that for 
A = 0.1, while PST can be qualitatively measured at t = tu, 
its subsequent periodicity is lost (not shown). This implies 
that in this case the gate should be put into effect at t = tu- 

Another potential source of gate error for the remote two- 
qubit gate is a timing error in the "extraction" of the qubits 
that have undergone the gate. Clearly from the parabolic ex- 
pansion of the concurrence around its maximum at t = tu, 
this error is second order in any timing error St, consistent 
with the PST timing error lllill . 
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Step 2: Injection at t=t /2 
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Step 3: Extraction and new injection at t=t 



FIG. 6: Diagram of steps 2 and 3 of the cluster ladder knitting proto- 
col. Solid lines indicate edges of the resulting cluster state 



V. DISTRIBUTED CLUSTER LADDER CONSTRUCTION 
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Having discussed the building block gate in detail, we now 
turn our attention to the repeated use of this operation to gen- 
erate a cluster ladder. A cluster state resource [ll6t between 
multiple qubits is realised by placing all qubits in state |+) 
and then applying a controlled-phase entangling gate between 
all pairs of qubits that are to be connected in the layout of the 
cluster state. The following protocol of injection and extrac- 
tion for a PST spin chain generates a cluster ladder: 

1. t = 0: Inject |+) at each end qubit (t = 1 and i = N) 
to a ground state chain (all zeros). This could be per- 
formed by a SWAP operation with register qubits adja- 
cent to each end site prepared in state |+). (Clearly the 
SWAP operation has to be fast on the timescale set by 
t M .) 

2. t = tnj/2: Inject |+) at each end qubit, which at this 
time will be disentangled from the rest of the chain (see 
Fig.©. 

3. t = tM- Extract the end qubits and inject |+) at each 
end qubit. 

4. Repeat last step as many times as desired, at time inter- 
vals of tj\i/2. 

5. Extract the last two qubits in the ladder when they reach 
the ends of the chain. 

It is assumed that each time a pair of qubits in the cluster 
ladder is extracted, each qubit is shuttled along a register at 
each end, with the next pair of | +) states moved into positions 
ready to be swapped in. This is illustrated in Fig. [6] 

Step 2 of the protocol, the second injection of |+) states, 
relies on the fact that the end qubits are initially empty. As 
especially for short chains, this cannot be guaranteed, we 
would then have some unwanted entanglement between the 
end qubits and the auxiliary qubits which carried the |+) 
states. If we assume that the injection is done by a SWAP 
operation, we can however refocus the system by measuring 
the auxiliary qubits: finding no excitations corresponds to a 
successful injection and also collapses the unwanted entangle- 
ment, whereas finding two excitations in the auxiliary qubits 
means that we have extracted all the excitations which the 
chain contained previously and are now back to the state of 



FIG. 7: Probability of successfully injecting at t = t m /2 vs. chain 
length iV. Data fitted as a function in 1/7V. 



the system as after step 1 of the protocol, allowing us to pro- 
ceed and attempt another set of injections a time later. 
In the event that we measure a single excitation in only one 
of the auxiliary qubits, the system cannot easily be recovered 
to a form which would allow us to continue with our protocol 
and needs to be re-initialised. 

The probability of successfully injecting a pair of | +) states 
at tnf/2 is however extremely high: Fig. Q shows the proba- 
bility of a successful injection with increasing N and we see 
that even for N = 9, the success probability is already next 
to 99%, with the failure probability (the deviation from unity) 
further decreasing as 1/N . These results allow us to be very 
confident about the efficiency of our proposed protocol. 

It must be noted that other injection methods, such as tun- 
nelling for example, might not allow for a refocussing mech- 
anism as described above. Injection methods which can result 
in one or both | +) states remaining in the auxiliary qubits will 
lead to retention of the unwanted entanglement between the 
chain and the auxiliary qubits as measurement of the auxil- 
iary qubits (which, in the case of a failed injection, are in a 
|+) state) does not lead to conclusive results. In this work we 
will therefore only consider injection of the SWAP type with 
refocussing taking place immediately after injection. 

An equivalent refocussing for the potentially imperfect ex- 
traction of |+) states at tM (step 3 of the protocol) is how- 
ever not possible. If the generated |+) states are extracted 
into empty storage qubits via a SWAP operation, while there 
is no unwanted entanglement between the chain and the stor- 
age qubits, there is no measurement we can perform to check 
whether the extracted states are indeed | +) states without de- 
stroying them. This might lead to some inaccuracies for very 
long cluster state ladders but results below suggest that this 
error does not limit practical use of the protocol. 

In order to now demonstrate our protocol, we consider the 
generation of a 'crossed' square cluster state, which simply 
involves injection of two pairs of qubits at t = and t = t^j /2 
with subsequent refocussing and their extraction at t = tM 
and t = 3tM/2 (see Fig. |SJ. First of all, we are dealing with a 
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chain of length N such that (N — 1) is a multiple of 4, so the 
gate in action in this case is, from Eq. ©, 

/l 00 \ 
G ,_ 1 

G - 10 0- (U) 
\ -1 / 

It is worth noting that G' is the product of a CZ gate, which 
leads to the desired entanglement, and a SWAP gate. In Fig. 
[H](a), this is illustrated by the labels on the lines representing 
excitations swapping as they cross whereas Fig. H](b) shows 
how the entanglement between the individual qubits it built 
up. Let us now look at the generation of a crossed square 
cluster state that this figure illustrates. 

To construct a crossed square cluster state, we inject exci- 
tations 1 and 2 at sites i = 1 and i = N respectively and then 
wait for tju/2 until we know that these sites are (nearly) com- 
pletely disentangled from the rest of the chain. We then inject 
excitations 3 and 4 at sites i = 1 and i = N respectively and 
refocus. The change in occupation probability of the individ- 
ual spins can be seen in Fig. [9ja) on the example of a 9-spin 
chain, for which the refocussing success probability is 0.9885. 

The occupation probability is defined as follows: let us as- 
sume that the set of k basis vectors forms the basis 
for our spin chain, such that any spin chain state \ip) can be 
written as = £\ =1 Cj|<fo), with £\ =1 |c 4 | 2 = 1. The ba- 
sis vectors can be represented as \<pi) = \j1j2 ■ ■ ■Jn}< with 
j = {0, 1}. Those basis vectors which contribute to the 
occupation of a site s are then \(f>i a ) = \j1j2 •■■ l a - - ■ j/v)> 
where 1 < s < N. Each \(j>i. s ) is weighed by its coeffi- 
cient Ci tS and so the occupation probability of a site s is given 

by 2»=i I c m| 2 - As a result, it is to be expected that the to- 
tal area of the histogram representing the occupation proba- 
bility is equal to the number of excitations in the system, so 
SsLi I c m| 2 = T for a specific excitation sector. 

In Fig. |9ja) we see how inaccuracies in manipulation of 
the spin chain lead to deviations from the ideal scenario. As 
the entropy of the depicted 9-spin chain does not reach zero at 
tu /2 (see Fig. |2}, the extremal spins 1 and N are not entirely 
decoupled from the rest of the chain. Despite refocussing af- 
ter injection, which guarantees an occupation probability of 
exactly 0.5 of spins 1 and N after injection at ijvf/2 (frame 
(a)), we see in frame (b) that the occupation probability of 
spins 1 and N before extraction at t = tm is bigger than 0.5. 
Consequently, in frame (c) the occupation probability of spins 
1 and N before extraction at 3£m/2 is smaller than 0.5. As 
we will see below, this will lead to a small loss in quality of 
the produced crossed square cluster state. 

After t = we have injected all the necessary exci- 

tations and wait for another 4m/2 until t = t^j before using 
a SWAP operation to extract excitations 3 and 4 at their in- 
jection sites, where the swapping qubits which are not part of 
the chain are presumed initially empty and now serve as stor- 
age. At this stage, the four excitations are already entangled, 
as can be seen in Fig. [8](b). Finally, at 3£m /2 we extract exci- 
tations 1 and 2 at their initial injection sites via another SWAP 
operation, thus completing the cluster state. Again, the auxil- 
iary qubits used to operate this SWAP are presumed initially 




FIG. 8: (a) Propagation of the four qubits involved in the genera- 
tion of a crossed square cluster state, (b) Schematic formation of the 
crossed square cluster state. 

unexcited and now serve as storage. The dynamics of these 
qubits is included in the simulation. Figure [9jc) illustrates 
that following this operation, the chain is virtually empty and 
thus ready for use again. We note that the residual occupa- 
tion in the chain is related to the imperfect disentanglement 
for a 9-spin chain of the end qubits at t = £j\//2 + nijvf , with 
n = 0, 1, 2, • • • . The end qubits disentangle virtually perfectly 
as longer chains are considered, see Fig. [2] 

Monitoring the quality of even just this four-qubit cluster 
state cannot be achieved via the EoF, which is only suitable as 
a measure of bipartite entanglement, but we instead consider 
the fidelity F against the ideal cluster state ipideai that we are 
hoping to achieve: 

F=\(i> tdea i\e- mt / h \ip ml )\ 2 7 (12) 

so perfect cluster state construction is achieved when F = 1. 
As in fact we are more interested in the true nature of the 
achieved state rather than simply in the amount of entangle- 
ment produced, this measure is very suitable for our pur- 
poses. Figure [TOl shows the evolution of this fidelity during 
and after implementation of the construction protocol for a 
crossed square cluster state, without the final read-out step at 
t = 3t m /2 so that the evolution of the achieved cluster state 
continues periodically beyond this time, achieving the same 
maximum fidelity at every 2t a/ after t = 'Hm /2. Notice the 
step at t = tjw/2, which corresponds to the injection of the 
second pair of excitations. This slight discontinuity is most 
clearly visible for N = 9, where we know from Fig. [2] that 
the dip in entropy is not as low as for longer chains. Nonethe- 
less, N = 9 achieves a fidelity of 0.9915 at t = 3t M / 2 and 
the two longer chains pictured both achieve unity at this time. 
This confirms that a 9-spin chain is long enough to demon- 
strate the effects we wish to highlight and so we will use this 
chain as our default device when considering a single value of 
N only. Refocussing for A = 13 and N = 17 is done with 
success probabilities of 0.9993 and 1 .000 respectively, show- 
ing that modest increases in N lead to both higher fidelity and 
higher refocussing success probability. As in Fig. |2] if no ex- 
citations are read out at t — 3tM/2, the evolution in Fig. [10] 
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FIG. 9: Changes in occupation probability of site i for a 9-spin chain 
during the generation of a crossed square cluster state. 
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FIG. 10: Fidelity of the ideal state vs. rescaled time t/tM for spin 
chains of length N = 9, 13, 17 used to build a crossed square cluster 
state. 




FIG. 11: Entropy S vs. rescaled time t/tM for spin chains of length 
N = 9, 13, 17, 25 used to build a crossed square cluster state. 



continues to be periodic over 2tM- 

We can confirm this result by also monitoring the entropy 
S, with the difference that instead of using the density matrix 
Pi.n of the two end spins, we use the density matrix of the 
four storage qubits, i.e. the two end qubits plus the auxiliary 
qubits. This is illustrated in Fig. QT| which clearly shows 
the entropy S dipping to zero (or virtually zero for N = 9). 
Again, the final read-out step at t = 3t^//2 is omitted, lead- 
ing to a continuous evolution of the constructed cluster state 
which displays dips in entropy every tu after the initial dip at 
t = 3<m/2. Analogous to the dip width decrease observed in 
Fig. |2] the variation in width of the dip at 1.5£m is approxi- 
mated by 1/y/N. 

S displays an additional dip at t = 0.75<a/, which is almost 
unnoticeable for short chains such as N = 9 but tends to zero 
for very long chains such as N = 25. This corresponds to 
a temporary decoupling of the end spins (sites 1 and N) as 
the excitations pass through each other (as indicated by the 
intersecting lines in Fig. |8). To illustrate this, we record the 
site occupation probabilities at t = 0.75tm in Fig. Q~2]for a 



short chain (N = 9 in panel (a)) and a long chain (N = 21 
in panel (b)). We see that for N = 9, where we see no dip 
at t = 0.75fM in Fig. QT| there is a high occupation prob- 
ability for the ends spins, whereas for N = 21, where we 
have a significant dip in entropy S at t = 0.75£a/, the occupa- 
tion probability of the end spins is very low, leaving the spins 
nearly decoupled from the rest of the chain. We also note that 
regardless of the chain length, the middle of the chain is empty 
at this time. 



VI. ERRORS IN THE CLUSTER KNITTING PROTOCOL 

Just like the two qubit gate in section II, the cluster knitting 
protocol we presented will be subject to various unpredictable 
errors and decoherences. Due to its higher complexity, the 
knitted cluster state is potentially more affected by defects 
similar to the ones discussed in section IV. Due to computa- 
tional restrictions, the more involved nature of the numerical 
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FIG. 12: Occupation probabilities of site i = 1, • • • JV at t = 0.75£a/ 
for a 9-spin chain (a) and a 21 -spin chain (b). 



simulations presented in this section also puts a limit on the 
lengths of spin chains we can investigate. The longest spin 
chain we consider here will therefore have 25 spins and not 
29 as in section IV. 

First of all, we re-consider the effect of random on-site en- 
ergies Ei. Again, each point in Fig. [T3land[T4lcorresponds to 
an average taken over 100 realisations. Figure [13] shows that 
the loss in fidelity still scales linearly with the number of spins 
in the chain, but the effect observed is much more detrimental 
than it was for the simple two-qubit gate (see Fig. [3j due to the 
increased number of excitations (now 4 vs. 2 in the two-qubit 
gate). Again, the decay in fidelity is not linear for increasing 
e I J m ax- While for very small values of e, fidelity is very well 
maintained even for very long chains, a perturbation of as lit- 
tle as 5% may already be very noticeable. Short chains up 
to TV = 17 achieve over 80% of the fidelity, whereas longer 
chains such as N = 25 suffer losses of over 30%. Similarly, 
when e = 0.1 acceptable fidelity is only achievable for very 
short chains such as N = 9, while chains of 17 spins and more 
become unsuitable for the proposed protocol as they achieve 
less than half of the desired fidelity. 

In Fig. Q3] we see the combined influences of non-uniform 
on-site energies E- t , weighted by e, and %' as given by Eq. (O 
added to the Hamiltonian @. Compared to Fig. [4] the effect 
of e has become extremely detrimental, leading to a loss in 
fidelity of over 60% for large values of e = 0.2. Even though 
a direct comparison of Fig. |41 and IT4lis not possible due to 
the different measures of state transfer quality, it is safe to say 
that the crossed square cluster state is much more affected by 
perturbations due to unwanted non-uniform on-site energies. 
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FIG. 13: Fidelity of the ideal state at t = 3t M /2 vs. N for three 
values of e, as labelled, where the chains are used to build a crossed 
square cluster state. Data points are averaged over 100 random on- 
site energy realisations. 
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FIG. 14: Fidelity of the ideal state at time t = 3tM /2 for a 9-spin 
chain used to build a crossed square cluster state vs 7 and e . Data 
points are averaged over 100 random on-site energy realisations. 



However, for smaller values of e up to a few percent, the fi- 
delity of the knitted crossed square cluster is still very good at 
about 90% of the ideal value. As the knitting protocol involves 
up to four excitations, the effect of Tl' is visible slightly more 
clearly but remains a very minor factor in the transfer quality, 
affecting the fidelity by a few percent only. The accuracy of 
Fig. |4] and [14] is the same, the perceived roughness of data 
points in Fig. [4] which is due to the randomisation of the in- 
fluence of £, does also appear to the same extent in Fig. [Pflbut 
is less visible due to the different scale. 

Finally, we also re-consider the influence of next-nearest 
neighbour interaction, as given previously by Eq. ([TOt . Fig- 
ure [15] shows that even for very small values of A below 5%, 
there is a noticeable loss in state fidelity. The decay scales 
now as N 2 for the shown values of A and N, resulting in fi- 
delities above 90% only for A = 0.01 for chains of lengths 
up to 25 spins. For chains of 9 spins a next-nearest neigh- 
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FIG. 15: Fidelity of the ideal state at time t = 3ijvr /2 vs. N for three 
values of A, as labelled, where the chains are used to build a crossed 
square cluster state. Lines are quadratic fits to guide the eye. 



bour coupling A of up to 5% still gives a fidelity of almost 
90% (not shown), while systems with 25 spins and A = 0.03 
lose the vast majority of their fidelity. For larger values of A, 
while the first fidelity peak at t = 'Mm /2 might still be of ac- 
ceptable value for short and medium length chains, it has to 
be noted that the periodicity is subsequently lost. Depending 
on the type of algorithm the obtained cluster state is intended 
for, this might pose a serious problem. For A = 0.05, chains 
of 21 spins or longer do not form fidelity peaks at t = Mm /2 
anymore, so that the protocol fails. This magnified detrimen- 
tal effect of next-nearest neighbour interaction is again due 
to the increase in number of excitations, a phenomenon that 
can also be confirmed in non-entangling spin chains subject 
to this perturbation (not shown). This underlines the necessity 
for strict control of longer range interactions, particularly in 
long chains. 

Another potential issue in the fabrication of distributed 
cluster ladders is that of untimely or non-synchronised injec- 
tion of excitations into the spin chain. The effect of the delay 
of a single qubit on a two-qubit gate has been discussed in Ref 
Ir29ll . but as the knitting of a distributed cluster state will in- 
volve four excitations most of the time, the possible delays are 
more involved. We will investigate the effect on the formation 
of a crossed square cluster state, the quality of which shall be 
measured via the state fidelity, and consider four separate de- 
lay scenarios. Referring to Fig. [8] the system is subjected to a 
delay St as follows: 

• (A): Excitations 3 and 4 are both injected at tm /2 + St 

• (B): Excitation 4 is injected at tnj/2 + St (while excita- 
tion 3 is on time) 

• (C): Excitation 2 is injected at St and excitation 4 is in- 
jected at + St (this corresponds to all injections 
on one side of the chain being delayed by the same 
amount) 
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FIG. 16: Effect of delay scenario A with St = O.Um on the fidelity 
of a 9-spin chain. 



1 




=a 0.75 - \- 

E 
< 

0.7 1 1 1 1 1 1 

0.02 0.04 0.06 0.08 0.1 



St/t M 

FIG. 17: Effect of delay scenarios A-D on the amplitude of the first 
fidelity peak of a cluster state on a 9-spin chain. 

• (D): Excitation 1 is injected at St and excitation 4 is 
injected at tj\//2 + St 

Refocussing as discussed in section V is done immedi- 
ately after each individual injection by measuring the auxil- 
iary qubits. Despite an initial perturbation to the generation of 
the crossed square cluster state, we see in Fig. [T6l(describing 
scenario A) that the evolution of the system continues to be es- 
sentially periodic if no read-out at 3f m /2 is undertaken. There 
is a clear kink in the plots of both entropy and state fidelity at 
£.a//2 + O.Um, where O.liM is the delay St of the injection 
of excitations 3 and 4. Consecutively, the fidelity peak occurs 
at Mm /2 + St and does not reach unity but a value of 0.9580, 
demonstrating the robustness of the system. 

Figure [T7] shows the amplitude of the first fidelity peak at 
3tM /2 + St for a range of St and all four scenarios considered. 
First of all we note that none of the delay scenarios lead to fi- 
delity of less than 90% for a delay time of 5% of the mirroring 
time tjvf, while delay up to 10% of 4m might lead to a loss of 
over 25% of fidelity. There is a clear discrepancy between 
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scenario A and the other three scenarios, with scenario A per- 
forming much better and suffering just over 10% fidelity loss 
for a delay of 10% of tu ■ Despite scenario B having fewer 
delayed excitations, it performs significantly worse. As such, 
it becomes clear that the system favours symmetrical input, 
also shown by scenario C with two excitations on the same 
end of the chain being delayed performing slightly better than 
scenario D, where excitations at opposite ends of the chain are 
delayed. 

VII. CONCLUSIONS 

In conclusion, we have presented a method to knit dis- 
tributed cluster states, using only a single spin chain set up 
for perfect state transfer and its natural dynamics. By closely 
observing the entropy of the chain, we have shown that the 
two end spins become decoupled from the rest of the chain 
at regular intervals, allowing us to inject further excitations 
into the chain without perturbing its existing excitation sub- 
spaces. The following natural chain dynamics lead the en- 
semble of excitations to entangle the end spins, which in turn 
decouple from the rest of the chain, allowing us to extract their 
states. As this routine of injections and extractions can be re- 
peated without theoretical restrictions, we are thus producing 
a knitted cluster state consisting of an even arbitrary number 
of spins. 

A further outlook of our technique is the possibility to knit 
other topological arrangements of cluster states by varying the 
injection and extraction timings and rates. As demonstrated in 
Fig. |6]and[8j each crossing of excitations leads to a new edge 
in the final produced state. Provided that the chain used to knit 
an arbitrary arrangement is long enough to allow for the nec- 
essary localisation of excitations, the number of entanglement 
bonds produced is therefore freely controllable. Examples of 
topologically useful entangled states that could be achieved 
in this way are presented in reference Q23J1 - We would how- 
ever expect more complicated structures which require very 
long chains, in particular those involving large numbers of ex- 
citations, to be more prone to the sources of decoherence we 
discussed. 

We also examined multiple potential causes of error when 
knitting states, both for the fundamental building block of our 



protocol, the two-qubit gate, as well as for the smallest knit- 
ted cluster state, the crossed square cluster state. Consistent 
with existing results on the influence of errors on state trans- 
fer in spin chains, the two-qubit gate is very robust against 
unwanted interactions between different excitations as well as 
perturbations in the on-site energies, maintaining high levels 
of over 90% of the desired entanglement even for large errors 
of 10% and very long spin chains. Longer range interactions 
on the other hand have potentially very detrimental effects, 
with large perturbations leading to a significant decay in en- 
tanglement as well as a loss of periodicity. Consequently, we 
observed analogous phenomena for the crossed square clus- 
ter state, amplified greatly by the increased number of excita- 
tions. Here, only non-uniformity in the on-site energies of up 
to 5% can be tolerated in order to ensure acceptable levels of 
over 70% of formed entanglement in long chains. Again, it is 
however the next-nearest neighbour interaction that deserves 
most attention in the fabrication process of a spin chain, as 
even small unwanted interactions of just a few percent can 
lead to significant fidelity losses, as well as loss of periodicity 
of the system. 

Additionally, we have considered the effect of the mistimed 
injection of one or more excitations on the formation of the de- 
sired crossed square cluster state and found that pairwise de- 
lay at the second injection time is much less detrimental than 
delay of a single excitation or delay of excitations at different 
injection times, but also that non-symmetric delay leads to a 
larger fidelity loss than delay restricted to the input on one end 
of the chain. Overall, the effect of delayed input as we con- 
sidered is however quite limited and does not lead to a loss of 
more than 10% of the perfect state fidelity for delays up to 5% 
of t M . 

Our studies have analysed and demonstrated a new protocol 
for the fabrication of distributed cluster states which requires 
no additional resources or mechanisms beyond the set-up of a 
spin chain for perfect state transfer and the associated SWAP 
operations used to inject and extract information. In the con- 
text of limited errors deviating less than 10% from the ideal 
set up, the knitting method proves to be a sound and stable 
method for generating distributed cluster states, making the 
utilisation of spin chains in quantum communication an ever 
more attractive prospect. 
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